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ABSTRACT 


^  The  modified  phase  space  theory  of  reaction  rates  and  Itonte 

Carlo  trajectory  calculate <kis  we  applied  to  the  problem  of  three- 
body  recombination  and  dissociation  of  diatomic  molecules.  Reac¬ 
tions  studied  include: 

N-  +  A  ^  U  +  N  +  Ar 
2  r 

HCf-  +  Ar  ^  H  +  C£  +  Ar 
+  H^H  +  H  +  H 

and  in  the  iai-nolecule  reactions 

F~  +  (At,  CO)  +  F  +  e"  +  (Ar,  CO) 

Br"  4  Ar  ^  Bi  +  e“  +  Ar 

On  the  basis  of  comparisons  .nade  between  the  phase  space  theory 
predictions  of  -.he  rate  ceofficients  and  e:rperimental  measurements 
ir.  is  concluded  tnet  the  pncse  space  theory  i*  capable  of  yielding 
good  a  priore  estimates  of  three-body  rate  coefficients  over  a  wide 
temperature  range.  ' 

Exploratory  experiments  were  carried  out  in  a  stratified- 
flow  wind  funnel  ir.  an  attempt  to  measure  the  growth  of  momentum- 
less  wake  in  a  stratified  flow. was  concluded  that  the  small 
wind  tunnel  was  not  a  sufficiently  refined  facility  for  making  the 
desired  measurements. 
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CHEMICAL  KINETICS 


Moote  Carlo  Trajectory  Calculations  of  the  Three-Body  Recombination 

(D* 

and  Dissociation  of  Dlatoiaic  Molecules 

The  modified  phase-space  theory  of  reaction  rates  was  ap¬ 
plied  to  the  problem  of  the  three-body  recombination  and  dissoci¬ 
ation  of  diatomic  molecules.  Tiie  results  illustrate  the  important 
influence  of  the  weak  attractive  minimum  in  the  third-body  inter¬ 
action  potential  and  the  effect  of  barrier  penetration  for  recom¬ 
bination  at  low  tenq)eratures.  The  system  ^2  +  Ar  was  used  as  a 
typical  illustrative  example,  and  Rood  agreement  was  obtained  be¬ 
tween  the  tiieoretical  predictions  and  the  experimental  measurements 
of  the  reaction  rate  coefficients  over  the  temyerature  range  200- 
12,000°K.  Tne  recrossing  factor  and  the  nonequilibrium  factor  were 
obtained  from  Monte  Carlo  trajectory  calculations  for  states  near 
the-  dissociation  limit.  Tlie  trajectories  were  sampled  v/ithin  the 
reaction  zone  wit’n  a  vjeight  proportional  to  the  equilibrium  reaction 
rate,  and  numerically  integrated  in  both  timewise  directions  to 
determine  the  complete  histories  of  the  collisions.  A  simple,  separ¬ 
able  fimction  for  the  equilibrium  transition  rate  R(e^,  £^)  from 
initial  energy  states  c.  to  final  energy  states  C^,  wliich  could 
;>e  ci<..racLerized  by  three  parnT-siters,  was  obtained  to  fit  the  numer- 
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Ical  (InCa  nrul  as  used  to  solve  tlie  steady-state  master  equation, 
ilstribut.lous  of  the  trajectories  v;ith  respect  to  energies  and  im¬ 
pact  parameters  are  presented,,  and  several  typical  trajectories 
are  shot/n  to  illustrate  the  important  features  of  the  collisions. 

The  contribution  to  the  over-all  reaction  rate  from  the  ’’complex 
meclianism  ’  was  also  obtained. 

(2) 

^!ontc^■  Carlo  Trajectory  Calculations  of  the  Dissociation  of  HC^  in  Ar 
The  modified  pahse-space  theory  of  reaction  rates  was  ap¬ 
plied  to  the  dissociation  of  HC2.  in  a  heat  bath  of  argon  atoms. 

Excellent  agreement  was  obtained  between  the  theoretical  predictions 
and  the  shock-tube  measurements  of  the  dissociation  rate  coefficient 
over  the  temoerature  range  2500-5000®K.  Tlie  recrosslng  correction 
factor  and  nonequilibrium  correction  factor  were  obtained  from  Monte 
Carlo  trajectory  calculations  for  states  near  the  dissociation  limit. 

Tlie  trajectories  were  sampled  within  the  reaction  zone,  with  a  weight 
proportional  to  the  equilibrium  reaction  rate,  and  numerically 
integrated  in  both  timewise  directions  to  determine  the  complete 
histories  of  the  collisions.  A  simple  separable  function  for  the 
equilibrium  transition  rate  R(Cj,,e^)  from  Initial  states  to  final 
states  was  obtained  to  fit  the  numerical  data  with  sufficient 
accuracy  and  was  used  to  solve  the  steady— state  master  equation. 
Important  features  of  collisions  of  highly  asymmetric  diatomic  molecules 

Preceding  Page  Blank 


are  discussed,  and  several  typical  reactive  trajectories  are  shovai 
to  illustrate  the  improtance  of  rotational  motion  in  such  collisions. 


Thermal  Dissociation  and  Recombination  of  Hydrogen  According  to  the 


Reactions:  112  +  H  H  +  H  +  H 


(3)  t 


Three-body  recombination  and  dissociation  rate  coefficients 
(k^  and  respectively)  for  the  reactions  H2  +  H  H  +  !1  +  H  have 
been  evaluated  for  temperatures  ranging  between  300  and  10,000®K 
on  the  basis  of  the  modified  phase-space  theory  of  reaction  rates 
and  :ionte  Carlo  trajectory  calculations.  The  semiempirical  Porter- 
Karplus  surface  for  was  used  In  the  calculations.  Good  agree¬ 
ment  was  obtained  between  the  theoretical  estimates  of  k,  and  the 

d 

bulk  of  the  dissociation  rate  measurements.  However,  in  contra¬ 
distinction  to  some  of  the  recent  experimental  investigations,  neitlser 
the  extremely  steep  temperature  decay  (ocT^ )  of  k  at  high  tempera¬ 
tures  nor  the  existence  of  a  pronounced  maximum  in  k^  at  moderate 
temperatures  are  supported  by  this  theoretical  work.  Possible  reasons 

for  these  differences  are  suggested.  A  simple  ratio >was  deduced 
/tj' 

which  relates  k^  ’  and  the  corresponding  rate  coefficients  of  the 
isotopes  of  hydrogen.  The  result  is  in  good  agreement  with  avail¬ 
able  experimental  measurements  and  enables  one  to  predict  the  latter 
rate  coefficients  without  additional  trajectory  calculations.  An 
example  is  given  to  illustrate  the  capability  of  the  present  Monte 
Carlo  method  in  obtaining  other  useful  information  about  atomic 
and  molecular  collision  processes. 


{■Research  initiated  under  ONR  Contract  No.  NOOOlA-67-0204-0040,  NR  061-141 
and  completed  under  NSF  Research  Grant  GP-32522X 
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Phase-Space  Theory  of  Electron  Detachirient  in  Slow  Atomic  Collislons  j 

A  classical  phase-space  theory  for  the  excitation  of  electrons 
in  slow  atomic  collisions  is  presented.  An  equation  which  governs 
the  energy  transfer  getween  the  electron  and  the  heavy  particles  is 
derived  for  slow  collisions  (adiabatic  approximation).  This  equation 
is  analogous  to  the  first  law  of  thermodynamics  for  an  open  system. 

The  theory  contains  only  one  adjustkble  parameter  and  g'ves  excellent 
correlation  for  existing  experimental  measurements  of  detachment 
rate  coefficients  of  F  +  (Ar,  CO)  and  Br  +  Ar.  Possible  mod¬ 
ifications  to  the  present  theory  are  suTunarized.  It  is  anticipated 
that  the  theory  should  be  applicable  to  other  electronic  rate  pro¬ 
cesses  such  ns  ionization  and  charge  trmsfer. 


'  Research  initiated  under  ONR  Contract  No.  N00014-67-020A-0040, 
NR  061-141  and  completed  mder  SSF  Research  Grant  GP-32522X 
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MOMHIUTL'ILIiSS  GROITOJ  IN  STRATIFIED  FLCR^ 

Experiirents  were  conducted  in  a  stratlfied-f low  -ind  tunnel 
to  deterinine  the  p.rowth  of  morcentumless  wakes  in  a  stratified  flow  . 
Zero  moioentum  wakes  vjere  fomed  behind  a  flat  disc  sup led  at  its 
center  with  an  axial  jet,  the  momentum  of  the  latter  bCiA  '  k justed 
to  balance  the  calculated  drag  of  the  disc.  Wake  grov/th  aii'  '•I'vi- 
tational  collapse  were  detected  visually  by  smoke  inj^ictioi  and 
by  sampling  with  a  hydrocarbon  analyzer  after  introducing  ethane 
into  the  jet.  Experiments  were  conducted  at  disc  Reynolds  numbers 
exceeding  200,  which  was  a  sufficient  condition  to  ensure  turbulent 
growtii. 

'.'aile  wal;e  growth  and  subsequent  collapse  were  visually  evi¬ 
dent  \;hen  smoke  was  introduced  into  the  wake,  the  iso-concentration 
lines  of  ethane  showed  no  marked  gravitational  collapse  correspond¬ 
ing  to  that  observed  with  the  smoke  tracer,  although  turbulent  wake 
growth  was  clearly  evident,  because  a  long  sampling  time  of  about 
a  minute  was  required  to  determine  each  concentration  reading,  it 
was  concluded  that  some  unsteadiness  in  vertical  flow  velocity  com¬ 
ponent  was  smearing  out  what  would  otherwise  be  steeper  vertical 
concentration  gradients,  causing  significant  differences  between 
the  near  instantaneous  smoke  trace  and  the  time-averaged  ethane 
concentrations. 

It  was  concluded  that  the  steady  flov?  technique  for  examining 
momentumless  wake  growth  and  collapse  in  a  stratified  flow  would 
encounter  significaric  difficulties  because  of  flow  unsteadiness 
whei)  operating  close  to  the  critical  Froude  number.  The  small  wind 
tunnel  used  in  these  experiments  was  not  a  refined  facility  and 
could  not  be  used  for  quantitative  measurements  of  the  type  desired. 
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